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Abstract 



i-G We use gauge/string duality to analytically evaluate correlation lengths of the renormalized 

^~-< field strength correlators in pure Yang-Mills theories at zero and finite temperature. 

Oh 

1 Introduction 
(N 
> 

(^ Understanding the infrared behavior of gauge theories from first principles is a longstanding prob- 

t^^ lem that offers perhaps the best hope of eventually understanding all the mysteries of quantum 

iJJ chromodynamics (QCD). It is well known that an important class of gauge invariant correla- 

• tion functions in the QCD vacuum is constructed from field strengths and Wilson lines. These 

(^ correlators do play an important role in several areas of QCD including models of stochastic 

<0 confinement of color, high energy scattering, heavy quarkonium systems, etcjj 
'T^ The simplest field strength correlator in four-dimensional Euclidean space is defined by 

> 
>^ D^,,p,(x) = (tr[G^,(0)C/p(0,x)Gp,(x)C/p(x,0)] ) . (1) 

d Here x^ are the Euclidean coordinates, the trace is over the fundamental representation, G^,/ 

is a field strength of the gauge field A^ and C/p(x,0) is a path-ordered Wilson line. The last is 
defined as C/p(xi, X2) = P eyi\>\ig Jq ds-^ Afj_{x{s))] , where s is a parameter of the path running 
from at X = xi to 1 at a; = X2 and 5 is a gauge coupling constantjj The paths are defined as 
straight lines. 

If one considers not a field strength but a field with different color quantum numbers, then 
the construction of correlators is changed. For instance, a field transforming in the fundamental 
representation of the gauge group will produce a two-point correlator 

^{x) = {q{0)Up{0,x)qix)), (2) 
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^The literature on the field strength correlators is vast. For a review, see, e.g., [T] and references therein. 

^In what follows we omit the indices when it is clear from the context. 



with q the field in the fundamental representation of SU{N) and q its conjugate. In general, 
there are many possible generalizations of the correlators by choosing different (curved) contours 
for the Wilson lines, or by inserting operators along the contours. 

It is believed that for large separation of the field strength operators the two-point correlator 
falls off exponentially [1] 

'^IJ.v^fyrix) ~ P^jy.pre"^ , (3) 

where Pfii^,pT is a kinematical factor, r = ^J{x^^, and A is a correlation length. This is the case 
also for ^, in fact 

*(x) ~ e"? , (4) 

with a correlation length ^. 

Until recently, the lattice formulation, still struggling with limitations and system errors, 
was the main computational tools to deal with non-weakly coupled gauge theories. The field 
strength correlators were also intensively studied (for a brief review, see p]). The situation 
changed drastically with the invention of the AdS/CFT correspondence |3j that resumed interest 
in another tool, string theory. The original duality was for conformal theories, but various 
perturbations (deformations) produce gauge/string duals with a mass gap, confinement, chiral 
symmetry breaking, etc jl]. 

In this note we continue a series of recent studies devoted to a search for an effective string 
description of pure gauge theories. In [5], the model was presented for computing the heavy 
quark potential at zero temperature. Subsequent comparison fH] with the available lattice data 
has made it clear that the model should be taken seriously. A non-trivial cross check for this 
model [7] , which checked the phenomenological value of the gluon condensate ^8j , was also car- 
ried out. Later, the model was extended to finite temperature. The results obtained for the 
spatial string tension |9j and the expectation value of the renormalized Polyakov loop [10] in 
the deconfined phase are remarkably consistent with the lattice, too. As is known, QCD is a 
very rich theory supposed to describe the whole spectrum of strong interaction phenomena. The 
question naturally arises: How well does such an effective string description capture other aspects 
of quenched QCD? Here, we address the issue of computing the field strength correlators in an 
analytical way as an important step toward answering this questionjj 

Before proceeding to the detailed analysis, let us set the basic framework. As in [5l [7| , we 
take the following ansatz for the five-dimensional background geometry 



BZ^ 



ds'^ = QNMdX^dX^^ = R^w {df + dx'^ + dz'^) , w{z) = ^ (5) 

to describe a pure gauge theory at zero temperature. The metric ([S]) is that of a deformed AdSs 
space, where 5 is a deformation parameter whose value can be fixed, for example, from the heavy 
quark potential. We take a constant dilaton and discard other background fields. 

When we go to the deconfined phase, we consider the five-dimensional geometry of [9l [TO] 



^To our knowledge, there have been no studies (numerical or analytical) of this issue from the viewpoint of 
AdS/CFT, or gauge/string duality, in the literature. 



ds' 



<5NMdX^dX^^ = R^w ( fdf + dx^ + ^dz"^ 



m = 1 - (^) 



(6) 



which represents a deformed Schwarzschild black hole in AdSs space. Here z-r is related to the 
Hawking temperature T = I/Ittz-y) whose dual description is nothing but the temperature of 
gauge theory. 

2 Calculating the Correlators at Zero Temperature 

To begin with, we need to find a recipe for computing the field strength correlators within 
the gauge/string duality. The strategy for finding it is as follows. In discussing the Wilson 
loops [IT], one first chooses a contor C on a four- manifold which is the boundary of a five- 
dimensional manifold. Next, one has to study fundamental strings on this manifold such that 
the string world-sheet has C as its boundary. The expectation value of the Wilson loop is 
schematically given by the world-sheet path integral {W{C)) = f DX e~^^ , where X denotes a 
set of world-sheet fields and Sw is a world-sheet action. In principle, the integral can be evaluated 
approximately in terms of minimal surfaces that obey the boundary conditions. The result is 
written as ( W{C) ) = ^^ w^e"'^", where Sn means a renormalized minimal area (in string units) 
whose weight is Wn- 
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Figure 1: Left: A minimal surface for a rectangular loop C indicated by thick lines along the (i,a;)-axes. 
It includes a lateral surface whose area is proportional to T and two identical end-surfaces with areas 
independent of T- Right: The end-surface in the x-z plane. It is bounded by a curved profile of a static 
string stretched between quark sources set at a; = 0, r and a straight Wilson line indicated by a thick line 
along the x-axis. Note that in the model we are considering a gravitational force prevents the string from 
getting deeper than z^ — 1/Vs into z direction j5j. 



In the case of a rectangular Wilson loop living on the boundary {z = 0) of five-dimensional 
space ([5]) a minimal surface corresponding to a static string cofiguration is sewn together from 
three surfaces, as shown in Figure 1 (left). For large T, we have |[5j 



{wicy 



-ET-2S 



(7) 



where E is a ground state energy of a quark-antiquark bound state and 5" is a renormalized area 
of the end-surface shown in Figure 1 (right). 



On the other hand, in the temporal gauge the expectation value of the Wilson loop can be 
evaluated as g 

{W{C))c^{tP\e-^'^\i;) = e-^^^\r), (8) 

where H is the gauge Hamiltonian in the temporal gauge and {ip) is the quark-antiquark bound 
state. It is given by \i>) = *(r) : g(r)C/p(r, O)g(O) : |0) such that (0| : q{r)Up{r,0)q{0) : |0) = 1. 
Since the quarks are non-dynamic, we need to normal order only the P-exponential function that 
results in non-zero expectation value for the normal ordered operator. Note that the exponential 
fall off Q is consistent with the fact that \ip) is normalizable. 
Then combining this with ([7]), we find 

^(r)~e-^, (9) 

where S" is a renormalized area of the surface shown in Figure 1 (right). Note that this formula 
is valid only for the ground state. 

It is worth noting that a similar representation has been suggested in [T3] for the expectation 
value of the Polyakov loop. In this case there is no need for the quark sources. The string is 
prevented from getting deep into z direction by a black hole geometry such that the maximum 
value of z is bounded by the black hole horizon. In the model with dynamic quarks it appeared 
in |14j . A difference here is that the Wilson loop goes along an internal direction. In general, it is 
natural to expect that in the presence of dynamic quarks the representation ([9]) breaks down at 
large separations due to string breaking. Physically, this means that quark bound states decay. 

To write a formal expression for the field strength correlator Q, let us think of the Wilson 
lines as forming a long, narrow rectangular loop in the x-t plane, as shown in Figure 1 (left) but 
with small T- The field strength operators are set at (0,0) and {r,0). The exponential fall off 
(pi) for large r is due to the Wilson loop expectation value. Subleading corrections come from 
quadratic fluctuations in the world-sheet path integral and from the field strength operators. 
They are expected to give a polynomial prefactor P in front of the exponential. Taking the limit 
T — )• 0, we have 

D^^,p,(r)~P^^,p,e-2^, (10) 

where S is the same renormalized area as in Q . Note that this formula only provides the leading 
exponent in the large r limit. 

In analyzing the formulas ([9| and (10), an interesting relation arises. Since S is proportional 



to r, we find that the correlation lengths ^ and A are related as 

A = ^«. (11) 

Given the background metric, we can calculate the renormalized area S by using the exact 
shape of the static string stretched between the heavy quark sources [5J. With the large r 
behavior in mind, it is technically suitable to add two pieces, shown in dashed lines in Figure 1 
(right), to the original surface. Their areas are each of subleading order in 1/r. As a result, the 
surface of interest becomes a rectangular in the x-z plane. 



'For example, see |12| and references therein. 



Now we are ready to use the Nambu-Goto action equipped with the background metric (pi) 

S = 2^ y dV ^det G^j^^d^X^dpX^' . (12) 

Next, we choose t^ = x and r^ = z. This yields 

5 = — r / dzw = / duw{u) , (13) 

TT Jo ^ ^c Jo 

where 3 = ^ and z^. = -p- We have also rescaled 2; as z = z^u. Since the integral is divergent 
at u = due to the z~'^ factor in the metric, we have to regularize it. We do so by imposing a 
cutoff e. 

To next-to-leading order in e, the integral is given by 



/ 



1 1 



duw = - + V^Erfi(l) - e + 0(e) , (14) 



e 



where Erfi(x) is the imaginary error function. We use the modified minimal subtraction scheme 
to deal with this integralrl So, we subtract the power divergence 1/e together with a constant 
c whose value must be specified from renormalization conditions. As a result, the renormalized 
area takes the forrrlj 

5= ■^Vs?-(c + ^/^Erfi(l)-e). (15) 

vr 

Combining this with (W| and ([9]), we get the correlation length ^ 

C=^(c + \/^Erfi(l)-e)'\ (16) 

It is now clear that T must go to zero faster than e. Indeed, in the limit T/e — t- the 



contribution of the lateral surface vanishes and, as a consequence, the relation ( 10 ) holds. With 
this fact understood, it is now straightforward to find the correlation length A 

A = -^(c + V^Erfi(l)-e)-\ (17) 

To define the model properly, we must specify the renormalization conditions. In doing so, 
we start with the correlation length ^. Our stringy construction suggests a natural condition 
^ = l/v^, where a is the string tension. In the model we are considering it is given by u = egs/vr 
[5]. Making an estimate requires some numerics. For SU{3) a value of g fixed from the heavy 
quark potential is g ss 0.62 |6]. This gives c ss 3.50. A typical value of s is s f« 0.45 GeV^ [6l[T5]. 
If so, a value of A is estimated to be 

X^0.20im. (18) 



^Note that the use of the modified minimal subtraction scheme in |5] allows one to adjust the value of a constant 
in the potential. Moreover, it also helps in finding a relation between the gauge/string duality result of [10] and 
the lattice data for the renormalized Polyakov loop. 

®Of course, the constant part being a number (y'7rErfi(l) — e) can be absorbed into c. However, it becomes 
temperature-dependent at finite temperature (as we will see in section 3), so we keep it. 



Is it a reasonable value? Although the lattice calculations of [16J claim a slightly bigger value 
A ~ 0.22 fm but with an error of order 0.03 fm, there is an obvious troublesome question. If the 
calculations are made in two different renormalization schemes, why are the results so similar? 
Unfortunately, we have no real resolution of this problem. 

On the other hand, choosing A = uiq , with mo ~ 3.64-^/(T the mass of the lightest glueball 
[T7] . gives c ^ 5.62. A simple algebra shows that in this case A w 0.13 fm. 

3 Calculating the Correlators at Finite Temperature 

According to [Ij, at finite temperature we must separately consider the electric and magnetic 
correlators [1]. So, we decompose G^i, into the electric and magnetic fields: Ei = Gq-i and 
Bi = ^EijkGjk- It is expected [IB] that at large separations the magnetic correlators show 
exponential fall off for any temperature, while the electric ones do so just below the critical 
temperature Tc. 

To study the magnetic correlator D)^ {x) = ( tr[i?i(0,x)f^p(x,0)i3j(0, 0)C/p(0,x)] ), we take 
the Wilson lines in the x-y plane and regard C as a long, narrow rectangle similar to that of 
Figure 1 (with t replaced by y). Thus, C is now a spatial Wilson loop. A crucial difference 
from temporal Wilson loops is that in the deconfined phase spatial ones obey an area law, with a 
spatial string tension crs{\ In other words, a (spatial) string stretched between two well-separated 
sources doesn't break. This allows us to use a formalism similar to that of section 2. So, the 
exponential fall off for large r is due to the spatial Wilson loop expectation value. The magnetic 
field operators and world-sheet fluctuations contribute to a polynomial prefactor P in front of 
the exponential. Taking the 7" — )> limit, we find 

D!^^(0~P..e-^^ (19) 

where 5* is a renormalized area of the surface shown in Figure 1 (right). 

Repeating the arguments of section 2, we add two pieces to the original surface to simplify 



further calculations. The Nambu-Gotto action, which we wrote before as (12), is now 

5 = ^ y dV ^det &^^,da,XNdpX^ , (20) 

where the background metric © is given by ([6]). Using the gauge t^ = x and r^ = z, it becomes 

S = ^rrdz^ = ^^ [\v^iv). (21) 

^ Jo V7 ^ztJo V7^' ^ ' 

Here we have rescaled z as z = z-^v. As in section 2, we regularize the integral over v by imposing 
a cutoff e. Note that in the deconfined phase the large distance physics is determined by the 
near horizon geometry of the deformed Schwarzschild black hole in AdSs space [9l [1^. So, the 
upper limit is Zt- 

To next-to-leading order in e, the integral is given by 



^Recall that as is not a physical string tension because it is not related to the properties of a physical potential. 



dv 



w 



^(4) 77 (_i. 1 1. iTl 
— n\^'^\ 4'4'2'4r2 



)+V^lW^Mhlhl$)+0{e), (22) 



r2r(3. 



where 1^2(0; b, c; a;) is the generahzed hypergeometric function and Tc is given by Tc = \/s/vr [9]. 
We use again the modified minimal subtraction scheme and look for a renormalized area. So, we 
have 



S = QTr 



^(4) TTH / 1.1 l.lZ?^ I A^ZkEill pa- 3 S.lT^N 
; "T yr 2 T^/3\1-^2U' 4' 2' 4 7^/' 



-1\ 1-^21^ 4, 4, 2' 4T^J 



(23) 



In the limit T/e — )• the contribution of the lateral surface is negligible. As a result, the 
correlation length is given by 



vr 



20V5r 



7r^-m-t'2 



ray 



'4' 4' 2' 4 T^ 



) + 



T^T 



T^T 



iF, 



2^4' 4' 2' 4T^ 



(24) 



To complete the picture, we must impose the renormalization conditions. Let us try A^ = 
l/(2-^/cjJ). The reason for this is that at zero temperature it is equivalent to ^ = l/\^ because 
the spatial string tension remains constant below T^ where it coincides with the (physical) string 
tension. In the model we are considering the temperature-dependence of the spatial string tension 
for T > Tc is given by [9] 



0".c 



(T^exp< 2 



1 



(25) 



where a is the string tension. Note that it is in a good agreement with the lattice data up to 
3-4Tc. 



A little experimentation with Mathematica soon shows that the expression (24), looking 



somewhat awkward, provides a very similar temperature dependence as the simpler expression 



(25). In Figure 2 we have plotted the results. 
Thus, we have 




(26) 



3.50. It is remarkable that for c ~ 3.22 the maximum discrepancy is 



with the previous value c : 
of order 0.2%. 

To study the electric field correlator D^fix) = {tv[Ei{0,x)Up{x,0)Ej{0,0)Up{0,x)] ), we 
take the Wilson lines in the x-t plane and regard them as forming a long, narrow rectangle 
shown in Figure 1. In the deconfined phase a static string stretched between two well-separated 
sources breaks. As a result, the dominant surface is given by a sum of two disconnected surfaces 
[13| whose renormalized areas vanish in the limit T — >• 0. This shows that in the deconfined 
phase there is no exponential fall off behavior of the electric field correlators at large distances. 
Such a conclusion looks attractive because it agrees with that of lattice simulations [19j . 
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Figure 2: Left: A comparison of different ^ curves for SU{N) gauge theory. The solid blue curve 



corresponds to (24 1 with c « 3.50 fixed at zero temperature. The dashed green curve represents the "best 



fit" with c ~ 3.22. In both cases, A = Xs\t=t^- The red dots denote w ^, with as given by (25 1. Right: 
for c « 3.22. 



A ratio ^/ 



4 Concluding Comments 

(i) We would like to emphasize again that the field strength correlators (II]) are scheme and path 

dependent. The former is due to a linear divergency of expectation values of Wilson loops. It 

is known as the perimeter law. The latter becomes clear if we let C be a circular loop of radius 

r/2. This leads to 5 ~ r^ rather than 5 ~ r at large r. Note that when a value of T is finite, 

the contribution of the lateral surface is no longer negligible. It matters for A in (JSJ). 

(ii) It is a common wisdom that magnetic and electric correlators are related to spatial and 

temporal Wilson loops, respectively. We have followed this philosophy in section 3 with our 

proposal for the two-point correlators. 

(iii) Because of scheme ambiguities, the correlation length may look not so good. It is therefore 

interesting that the difference 



A-i(ri) xjHt2 



2g 



Ti T2 

is finite. Performing the integral gives 



dv 

y2 



{i-v^y 



exp 



rp2 



exp 



{ 



rp2' 



(27) 



\-\Ti) \-\T2 



Ti 



T2 



2^Q 



+ 



r(|) 

n\) 

r(|) 
r(i) 



1F2 



rp2 



1. 1 

'4' 4' 



2 ' 4 r| 



) 1-^2;^ 41 4> 2' 4T2y' 



4 ' 4 ' 2 ' 4 Tf - 



M\ 



4)4)2' 



3. l^c^ 
477 /* 



^2 



M\ 



_ 3 3. 1^2 

4' 4' 2' 4 7^ 



(28) 



Thus, our model predicts the scheme-independent relation between the correlation lengths of 
magnetic operators at different temperatures. It will be interesting to see whether it will match 
or close to match lattice simulations. 

(iv) If we consider not a field strength but an arbitrary local operator in the adjoint representation 
of SU{N), then it follows from our proposal that the correlator shows the exponential fall off 



(C(0)C/p(0,x)C'(x)C/p(x,0)) ~e-A for r ^ oo , (29) 

with a universal correlation length A. 
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